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19.  ABSTBACT  (Continued) 


after  Initial  Injection,  and  that  microscale  and  — processes  do  aot  break 
up  or  remove  the  aoot  layer The  validity  of  these  assumptions  la  discussed. 

Global  climate  aodelsj  (general  circulation  ends  Is  and  energy  balance  models) 
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a  few  days  after  Initial  Injection,  and  that  slcroscale  and  aasoscale  processes  do 
not  break  up  or  remove  the  soot  layer. j  the  validity  of  these  assumptions  is  dis¬ 
cussed  <^Xheuristlc~anaiysis~taitlt*tes  that  a  horizontally  uniform  and  optically 
chick  soot  layer  Is  thermally  unstable  and  will  break  up  Into  long  vertical  filaments 
This  process  aay  significantly  change  Che  seen  optical  properties  of  the  cloud  layer 
and  thus  alter  the  synoptic  and  planetary  scale  effects  of  the  soot  Injection. 
Suggestions  are  presented  for  further  Investigations  which  can  reduce  the  lar^~uncer 
talntles  associated  with  the  nuclear  winter  problem. 
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Tbrm  types  of  long- ter*  environments!  effects  9f  nudes?  weapons  bent 
ben  widely  debated — fallout.  ozcne  layer  depletion,  end  cllsats  Modifi¬ 
cation.  The  Issue  of  biological  daaage  from  radioactive  fallout  was 
considered  in  a  1975  study  by  the  national  today  of  Sciences  (MAS) . 

Their  report  (Reference  1)  concludes  that,  although  local  fallout  would 
constitute  an  luaedlats  danger  to  the  unsheltered  population  directly 
dcaeerind  fro*  a  target  eras,  long-term  and/or  wldaapraad  effects  would  be 
minor  bscauss  of  atmospheric  dilution  of  'the  radioactive  oonponant,  leading 
to  greatly  reduced  explosure  levels,  Sswever.  a  ana  recent  report  (Refer¬ 
ence  2)  concludes  that  the  global  fallout  problem  ay  be  aore  serious  than 
estimated  by  the  NAS.  though  typical  doses  would  still  ba  well  below  the 
lethal  range.  ■  ■ 


The  na jar  oeneem  of  ths  1975  NAS  report  was  the  affect  of  bemb- 
producad  oxides  of  nitrogen  an  stratospheric  aeon e  (Reference-  3) .  In- 
erseasd  .levels  of  ultraviolet  radiation  trarsmlttad  by  a  depleted  aeons 
layer  could  be  a  long-term  and  globel  danger  to  ths  survlvlng  populatlcn. 
However,  the  analyeis  by  the  MAS  asouned  Individual  weapon  yields  In  the  1 
to  20  KT  range.  The  subsequent  shift  in  the  O.S.  and  Soviet  arsenals  to 
smaller  and  nore  accurate  imapona  has  reduced  this  concern  .  since  plumes 
from  tropospheric  explosions  with  yields  below- 1  MT  are  unlikely  to  rise  to 
the  height  of  the  ozone  layer  (Reference  4) . 


The  third  type  of  effect— climate  modification  has  recently  received 
a  great  deal  of  attention.  In  1982.  Crutzan  and  Bizis  (Reference  5)  hypo¬ 
thesized  that  widespread  forest  fires  ignited  by  naclaar  weapons  would 
crests  a  thick  smote  layer..  This  would  redoes  ths  amount  of  sunlight 
teaching  ths  ground.  In  1983.  Tlxrco,  art  ad. .  (Reference  6)  argued  that 
urban  fires  would  create  an  even  aore  opaque  layer  of  black  soot,  vtfiich 
would  remain  In  the  atmosphere  for  several  weeks  or  longer.  During  this 
tine,  the  soot  could  spread  over  ouch  of  the  Northern  Hemisphere.  By 
absorbing  incoming  sunlight  well  above  ground  level,  thus  short-circuiting 
ths  normal  greenhouse  effect,  ths  soot  layer  would  lead  to  subfreezing 
surface  temperatures,  even  Is  summer. 


A  1986  study  by  ths  NAS  (Reference  7)  concludes  that,  although  the 
uncertainties-  in  the  TUrco.  et.  al. ,  scenario  are  presently  too  large  to 
allow  a  quantitative  prediction,  the  possibility  of  a  severe  surface  tem¬ 
perature  reduction  following  a  major  nuclear  conflict  should  be  seriously 


considered.  Such  a  climate  effect  wculd  damage  crops  and  lead  to  Increased 
casualties  among  survivors  of  Immediate  nuclear  effects  (Reference  8) .  The 
altered  atmospheric  conditions  could  also  have  synergistic  effects  on 
fallout  and  ozone  depletion  (Reference  2).  In  the  case  of  a  protracted 
conflict,  reduced  atmospheric  transmission  and  increased  pollutant  loading 
would  seriously  (apart  on  the  effectiveness  of  surveillance  assets  and 
air-breathing  weapons  systems.  For  the  above  reasons,  it  Is  desirable  to 
reduce  the  uncertainties  In  the  climate  predictions  to  a  level  vtoere  the 
public  policy  and  military  Implications  can  be  clearly  understood. 

Three  types  of  uncertainties  are  associated  with  the  climatic  pre¬ 
dictions:  the  nuclear  exchange  scenario  (number  and  size  distribution  of 
warheads,  target  selection,  ground  bursts  vs.  air  bursts,  and  ambient 
meteorological  conditions):  smoke  production  and  Injection  (amcunt  of 
smoke,  size  distribution  and  optical  properties  of  smoke  particles,  height 
of  the  smoke  plume,  and  aging  of  the  smoke) ;  and  atmospheric  response 
(radiative  effects,  circulation  changes,  smoke  transport,  convection  and 
precipitation,  and  cloud  interactions;  The  first  type  of  uncertainty 
requires  that  climatic  response  be  evaluated  as  a  function  of  the  type  of 
conflict  and  of  the  time  of  year  at  viiich  the  conflict  occurs.  The  effec¬ 
tiveness  ol  any  defensive  systems  must  also  be  considered.  The  second  type 
of  uncertainty  can  be  addressed  by  experimental  and  theoretical  Investi¬ 
gations  of  fire,  smoke,  and  plume  properties.  The  last  type  of  uncertainty 
must  ultimately  by  addressed  by  cenputer  simulation,  since  neither  natural 
nor  controlled  climate  perturbati/jns  occur  on  the  estimated  scale  of  the 
nuclear  war  effects.  In  the  absence  of  an  ability  to  observaticnally 
verify  the  simulations,  a  residual  uncertainty  regarding  the  predicted 
atmospheric  response  will  always  remain.  Some  specific  issues  can, 
however,  be  addressed  observational ly  by  studies  of  transport  and  aging  of 
smoke  frcar.  natural  and  man-made  fires. 

This  report  is  concerned  with  the  third  type  of  uncertainty,  and 
concentrates  on  modeling  problems.  Section  2  describes  simulations  with  a 
simple  type  of  climate  model.  The  simulations  were  carried  out  to  demon¬ 
strate  the  feasibility  of  using  this  model  to  test  atmospheric  response  to 
soot  loadings.  The  model  is  used  to  estimate  sensitivity  of  the  predicted 
cooling  to:  (1)  the  season  at  which  the  soot  is  injected.  (2)  the  soot 
scavenging  rates.  (3)  the  geographical  distribution  of  the  injected  soot, 
and  (4)  the  Infrared  properties  of  the  soot  particles.  Section  3  summa¬ 
rizes  the  modeling  results  and  provides  analyses  of  some  processes  not 
accounted  for  in  present  atmospheric  simulation  codes.  Of  particular 
Importance,  it  appears  that  the  uniform  soot  layer  assumed  In  these  models 
is  ttermally  unstable  and  will  tend  to  break  up.  The  Appendix  provides 
derivations  of  sane  mathematic  properties  of  the  climate  model. 
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SXKOXIGNS 


Climate  modeling  has  progressed  rapidly  oner  the  past  decade  due,  in 
part,  to  the  application  to  eliaata  prohlaaa  of  general  circulation  models 
(GCM) .  A  GCM  vmea  a  discrete  grid  aeeh  (or  the  equivalent  harmonic  expan¬ 
sions)  to  solve  the  primitive  equations  governing  dynamical  circulation, 
radiative  tranafer,  and  some  form  of  a  hydrological  cycle.  A  typical  GCM 
grid  has  horizontal  resolution  of  *  900  lot  and  vertical  resolution  on  the 
order  of  1  to  10  loa. 


Many  atmospheric  processes  occur  on  spatial  scales  smaller  than  the 
GCM  resolution.  Examples  of  subgrid  processes  include:  convection, 
precipitation,  cloud  interactions ,  frontal  zone  circulations,  and  boundary 
layer  effects.  These  processes  are  rapreeented  in  the  GCMs  in  highly  para¬ 
meterized  forma.  The  representations  contain  adjustable  coefficients  vhich 
are  "tuned"  to  produce  a  good  simulation  of  present  climatic  conditions. 


Climate  models  can  be  tasted  by  comparing  their  predictions  to  recent 
minor  climate  perturbations,  such  as  caused  by  volcanic  dust  veils,  and  to 
paleoellmates  reconstructed  from  proxy  data,  only  the  latter  type  of  test 
involves  departures  from  the  present  climate  of  a  magnitude  similar  to  the 
changes  postulated  to  occur  in  the  aftermath  of  a  imrlsar  war.  For  ex¬ 
ample,  a  few  interesting  ice  age  experiments  have  been  carried  out.,  with 
drastic  departures  from  the  present  climate  (References  9  and-  10) .  How¬ 
ever.  Interpreting  the  results  of  such  experiments  is  difficult,  due 
both  to  ambiguities  of  the  paleocliaate  data  base  and  to  the  inherrent 
complexity  of  GCM  diagnostics.  The  model  climate  generated  by  a  GJM  is 
often  ae  complex  and  inscrutable  am  the  real  climate. 


In  parallel  with  GCM  construction,  a  number  at  simpler  models  have 
been  developed  for  sensitivity  experiments.  It  Is  hoped  that  these  models 
will  produce  climates  which  respond  tn  external  influences  in  much  the  same 
war/  as  does  the  real  climate.  If  so.  their  simplicity  can  lead  to  a  better 
understanding  of  climate  response,  as  well  as  allowing  efficient  evaluation 
of  a  range  of  scenarios.  In  many  ways,  these  simpler  models  offer  a  good 
match  to  the  crudeness  of  the  present  understanding  of  the  type  of  pertur¬ 
bation  which  would  result  from  a  large-scale  nuclear  war. 


One  example  of  such  a  simple  model  is  a  radiative-convective  model 
(RCM) ,  such  as  used  in  the  study  by  Turco,  at  al.,  (Reference  6).  This 
type  of  model  has  excellent  vertical  resolution  and  handles  radiative 
transfer  in  considerable  detail  tut  has  no  horizontal  structure. 
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A  different  type  of  climate  model  is  an  energy  balance  model  (EBM) , 
popularized  by  Budyto  (Reference  11)  and  Sellers  (Reference  12).  and 
recently  reviewed  by  North,  et  al. ,  (Reference  13).  An  HMB  uses  very'  sim¬ 
ple  representations  of  energy  fluxes  to  examine  how  these  fluxes  interact 
to  determine  the  surface  temperature.  Most  EBMs  include  seme  horizontal 
resolution,  so  the  effects  of  geography  (land-sea  fractions)  are  included 
in  the  models.  In  a  sense,  EBMs  complement  RCMs,  which  directly  calculate 
energy  fluxes  but  have  as  geography. 

With  cne  exception  (Reference  14) ,  EBMs  have  not  been  applied  to  the 
Turco.  et.  al..  hypothesis.  The  remainder  of  this  section  describes  the 
application  of  a  continent-resolving  EEM  to  estimating  the  climate  impact 
of  a  soot  layer  of  the  type  postulated  to  follow  a  major  nuclear  conflict. 


sracxricancH3  op  the  model 

Me  use  a  two-dimensional,  one-level,  time-dependent  EBM.  of  the  type 
discussed  by  North,  et  al..  (Reference  IS).  This  model  wes  originally 
developed  for  studying  tne  initiation  of  ice  ages  following  changes  in  the 
seasonal  distribution  of  sunlight  (due  to  changes  in  the  Earth's  orbit). 

The  model  predictions  for  this  case  are  in  excellent  agreement  with  the 
regular  cycles  of  ice  ages  recently  discovered  in  the  paieoclimatic  data 
(Reference  16) .  The  model  represents  the  large  difference  in  thermal 
response  times  of  the  land  and  oceans  by  differing  heat  capacities.  The 
seasonal  cycle  (but  rot  the  diurnal  cycle)  is  incorporated,  so  the  seasonal 
dependence  of  nuclear  war  effects  can  be  studied.  Finally,  the  model  is 
simple  trough  that  its  principal  mathematical  properties  can  be  directly 
evaluated,  and  the  model  can  be  economically  exercised  for  a  variety  of 
cases.  The  price  pale  for  this  convenience  is  that  certain  atmospheric 
processes  (radiative  transfer  and  horizontal  energy  advectian)  are  treated 
in  highly  idealized  forms.  These  simplifications  can  beat  be  justified  by 
the  fact  that  the  model  does  predict  a  climate  v*iich  is  a  good  facsimile  to 
the  real  climate. 

The  basic  energy  balance  equation  is 

^  3T(r.t)  -  .  - 

C(f)  -  -  Q  a(x)  S(x,t)  -  [A  +  B  T(r,t) ]  +  DV2  T(r.t),  (1) 

at 

where  , 

t  -  time,  in  years  with  t  ■  0  at  winter  solstice, 
r" *  geographical  position  (latitude  and  longitude), 
x  a  sine  of  latitude, 

C(r)  »  heat  capacity  per  unit  area. 
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T(r»;)  *  surface  temperature  (dependent  variable) . 

0  ■  solar  ccnstaot/4. 

a(x)  -  fraction  of  solar  energy  absorbed  by  the  atmosphere- ocean 
system, 

S(x.t)  •  relative  solcr  radiation  distribution  vs.  latitude  and 
w®soo# 

A.B  *  empirical  infrared  radiation  coefficients,  and 

D  ■  diffusion  coefficient  f^r  horizontal  energy  transport. 

In  simple  terms,  Equation  (1)  states  the  energy  conservation  rule  for 
a  vertically  averaged  column  of  ataoephe>~c; 


net  rate 
of  heating 


absorbed 

sunlight 


emitted 

infrared 


horizontally 
transported  heat 


The  numerical  coefficients  in  this  model  are  derived  from  physical 
arguments  (C),  from  measurements  (Q,  a.  A,  and  a),  ae  by  "tuning"  the  model 
to  reproduce  tne  observed  climatology'  (D) .  Table  1  gives  the  values  used 
in  the  model  calculations  described  in  this  report:  justification  for  these 
values  is  discussed  belcw. 


The  effective  heat  capacity,  C.  depends  an  the  mass  of  air  and  water 
which  is  closely  coupled  to  surface  conditions.  Over  land,  we  use  half  the 
mess  of  a  colunn  of  atmosphere  (i.e. .  the  atmosphere  below  the  500-mb  level 
or  roughly  below  3  tan  altitude) .  Over  oceans,  we  add  in  the  heat  capacity 
of  the  water  in  the  wiaxi-driven  mixed  layer  of  the  oceen,  taken  to  be  uni¬ 
formly  73  meters  in  depth. 

The  solar  constant  (radiant  energy  flux  per  unit  area  at  the  top  of 
the  sunlit  atmosphere)  is  a  directly  measured  quantity  (Reference  17) :  the 
diurnal  average  is  obtained  by  dividing  the  solar  constant  by  four.  The 
fraction  of  the  solar  energy  absorbed  is  taken  to  be  a  quadratic  function 
of  latitude: 

a(x)  -  ao  +  a^x2.  (2) 

with  coefficients,  derived  from  zanally  averaged  measurements  of  direct  and 
reflected  sunlight,  obtained  fron  satellite  data  (Reference  18).  The 
Infrared  coefficients.  A  and  B.  are  likewise  derived  by  correlating  sate¬ 
llite  measurements  of  out-going  infrared  fluxes  with  surface  temperature 
(Reference  19) . 
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Tabic  1.  Model  coefficients 


Quantity 


Heat  Capacity 

C 

land:  5.07  x  10®  J/m2/°C 

oceans:  3.04  x  103  J/m2/°C 

Average  Insolation 

.  Q 

340  W/m2 

Co-Albedo 

a 

0.73 

a 

-0.18 

Infrared  Coefficients 

A. 

203.3  W/m2 

B 

2.094  W/m 2/°C 

Di “fusion  Coeffir lent 

D 

0.35  W/°b 

Seasonal  Cycle 

•  S1 

-0.796 

S2. 

-0.477 

So-, 

+0.147 

The  diffusion  coefficient  Is  not  directly  related  to  •  Measurable 
quantity.  Therefore,  the  value  of  this  parameter  was  adjusted  to  obtain 
agreement  of  the  predicted  and  observed  pole  to  equator  temperature  differ¬ 
ence  averaged  over  the  sewn-ial  cycle.  This  is  the  only  "ttned"  parameter 
In  the  model. 

War  simulating  the  normal  climate,  the  model  Is  driven  through  the 
seasonal  cycle  by  the  time-varying  distribution  of  sunli^it,  given  by 

S(x.t>  -  1  +  SjCOS(2*t)  P1(x)  +  (S2  +  S22  coa(4*t)>  P2(x).  (3) 

where  P.  and  P2  are  the  Legendre  polynomials.  The  coefficients  In  this 
expansion  were  adjusted  to  provide  a  "best  fit”  to  the  sore  complicated 
exact  expression  for  the  swasmel  cycle  (Reference  19). 

The  primary  mathematic  1  properties  of  the  model  concern  the  e-foldlng 
time  far  surface  cooling,  then  the  solar  heating  Is  removed,  and  the  con¬ 
duction  length  *  ale  for  propagation  of  a  local  temperature  perturbation. 
Derivations  of  these  properties  are  provided  In  the  Appendix.  Hie  results, 
with  the  coefficients  In  Table  l.  are 

Cooling  Time:  x  “  C/B 

■  28  days  for  continents 

■  4.6  years  for  oceans 

Conduction  Length:  &  ■  D/B 

-  4,500  km. 

In  terms  of  understanding  the  results  of  the  simulation,  we  note  that  the 
very  different  (by  a  factor  of  60)  cooling  times  for  land  vs.  ocean  imply 
that  continental  areas  will  cool  much  faster  than  oceanic  regions.  The 
long  conduction  length  Implies  that  the  stabilizing.  Influence  of  the  oceans 
will  propagate  well  into  the  continents.  This  Is  In  agreement  with  GCM 
results  (Reference  20) .  which  predict  that  the  most  severe  temperature 
effects  are  restricted  to  the  deepest  continental  interiors. 


MODEL  9QL0TICW  AMD  VALIDATION 


To  solve  the  energy  balance  equation,  the  spatial  dimensions  are  first 
expanded  Into  spherical  harmonics: 


L 

E 


l 

Z 


Awd  m*0 


£ra 


(c) 


Vr> 


T(r.t) 


(4) 


and 


C(r) 


2 L 

I 

tmQ 


* 

i  0^Vr' 


/ 

/ 


(5) 


where 


Y;ja(r)  -  P^fx)  ccs(wt>)  (6) 

and  $  is  longitude,  the  coefficients  C*.  are  evaluated  from  the  specified 
land-sea  distribution.  For  the  results  alscuawed  in  this  report,  we  used 
L  «  11,  implying  a  15°  x  15°  spatial  resolution  (abort  1,000  Ion  at  mid- 
latitudes)  . 

Substitution  of  Equations  (2)  through  (5)  into  Equation  (1)  gives  a 
system  of  (L  +  l)(L  +  2)/2  ordinary  differential  equations  for  the  mode 
amplitudes  T  (t).  These  equations  are  integrated  frcm  prescribed  Initial 
values  using™  fourth-order  Hamming  predictor-corrector  method.  Tempera¬ 
ture  maps  are  recooposed  and  plotted  at  desired  Intervals.  To  give  sane 
feeling  for  the  computing  requirements ,  simulation  of  1  year  of  climate, 
with  temperature  maps  at  10-day  Intervals,  requires  “  30  minutes  of  CFO 
time  on  a  VAX  11/750  ccnputer. 

The  ability  of  the  model  to  simulate  atmospheric  processes  was  tested 
in  two  ways.  First,  the  predicted  seasonal  cycle  of  an  unperturbed  (con¬ 
trol)  model  was  compared  to  the  observed  cycle.  This,  teat  is  not  redundant 
with  the  tuning  process  since  the  model  was  tuned  to  match  the  annual  mean 
climate,  not  the  seasonal  variation.  The  second  type  of  test  involved 
comparisons  of  EBM  predictions  for  perturbed  (soot)  models  with  results 
obtained  from  GCMs,  using  identical  soot  parameters.  This  is  not  a  direct 
validation  since  the  GCM  results  could  be  incorrect.  However,  we  do  not 
claim  to  be  able  to  make  predictions  v*iich  are  superior  to  the  GCM  results. 
The  advantage  of  the  EEM  lies  in  its  simplicity  and  its  ability  to  economi¬ 
cally  evaluate  many  scenarios .  with  integrations  covering  an  extended 
period  of  time. 

Figure  1  compares  the  amplitude  of  the  model  seasonal  cycle  with  the 
observed  climatology,  as  reported  by  Crutcher  and  Meserve  (Reference  21) 
and  Taljaard.  et  al. ,  (Reference  22).  The  observational  data  have  been 
smoothed  with  an  L  *  11  barmcnic  filter  (i.e..  the  data  were  analyzed  as  in 
Equation  ( 4 ) ,  and  the  map  in  the  bottom  panel  of  Figure  1  was  reconstructed 
from  the  derived  harmonic  coefficients) ;  this  procedure  sets  the  spatial 
resolution  of  the  data  equal 'to  that  of  the  model.  Both  panels  clearly 
show  the  Influence  of  land  masses,  with  reduced  heat  capacity  and  larger 
seasonal  amplitudes,  as  compared  to  the  oceans.  A  closer  examination 
Indicates  that  the  model  amplitude  tends  to  be  a  bit  small  over  North 
America  and,  perhaps,  sanevPat  too  large  over  Eurasia. 
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Figure  2  onnpnrcs  the  iinennnil  phase  lags,  defined  as  the  tin  between 
merlin  hi  solar  heating  and  temperature.  Results  within  +15  degrees 

of  the  equator  are  not  ahown.  because  the  amplitude  of  the  cycle  is  small 
In  this  region:  so  the  phase  lags  are  poorly  defined.  Both  the  model  and 
the  observations  shew  phase  lags  of  the  order  of  30  days  In  continental 
areas,  consistent  with  the  cooling  time  scale  discussed  earlier.  Over 
oceans,  the  phase  lag  in  the  observations  is  roughly  75  days.  The 

aodel  tends  to  produce  a  nearly  uniform  phase  lag  of  90  deys  (one  quarter 
cycle)  over  oceans.  This  corresponds  to  the  lag  of  a  linear  system  In  the 
limit  of  Infinite  heat  capacity.  Considering  the  very  email  amplitude  of 
the  oceanic  cycle  (both  in  the  observations  and  in  the  model) ,  we 

do  not  view  the  discrepancy  in  the  phase  lags  as  a  serious  problem. 

The  second  pert  of  the  validation  procedure  involved  repeating  soot 
experiments  **iich  had  been  run  cn  GCMs.  The  GCM  results  were  taken,  in 
graphical  farm,  from  published  reports.  Because  we  did  not  have  the  GCM 
results  available  in  digital  format,  it  wee  not  possible  to  perform  the 
spatial  filtering  and  replotting  v*iich  would  be  desirable  for  direct  com¬ 
parison  with  the  2EM  results.  A  second  aspect  of  these  comparisons  con¬ 
cerns  the  mode  in  which  the  JEM  was  run.  For  ail  soot  experiments,  the  EBM 
used  a  land-sea  distribution  synmetrized  about  the  prime  meridian.  Thus,  a 
Western  Hemisphere  simulation  actually  refers  to  a  globe  with  the  western 
continents  repeated  at  180  degree  intervals  of  longitude,,  and  similarly  for 
Eastern  Hemisphere  simulations Experiments  showed  that  this  type  of  east- 
west  symmetric  model  produced  essentially  the  same  results  as  an  asynmetrlc 
model  but  with  much  less  computing  time  since  only  the  even-m  modes  of  the 
spherical  harmonics  are  used.  Hie  effect  an  the  results-  is  primarily  cos¬ 
metic  in  that  the  eastern  and  western  maps  do  not  join  smoothly  at  the 
bexirdary. 

Covey,  et  al..  (Reference  20)  used  the  National  Center  for  Atmospheric 
Research  (WEAR)  Ccranamity  Climate  Model  (CCM)  to  simulate  the  effect  of  a 
fixed  soot  layer.  The  CCM  solves  the  general  circulation  equations  using 
spherical  harmonics  with  an  L  »  15  cut-off.  Vertical  resolution  is  pro¬ 
vided  by  a  finite  difference  scheme  with  9  levels.  Sea  surface  tempera¬ 
tures  are  fixed,  so  no  ocean  response  is  permitted.  The  swot  parameters 
assumed  by  Covey,  et  al.,  are: 

(1)  purely  absorbing  at  visible  wavelengths, 

(2)  optical  depth  *  3.0  at  visible  wavelengths, 

(3)  latitude  of  soot  bend  *  30°  to  70°N,  and 

(4)  no  effect  an  Infrared  radiation. 

Figure  3  shows  the  surface  temperature  changes  (with  respect  to 
control  runs)  for  simulations  starting  an  March  22  (top  panel)  and  June  30 
(bottom  panel) .  In  both  cases,  results  are  averaged  over  a  5-day  period  to 
smooth  the  meteorological  "noise"  of  the  model.  The  contours  shew  inter¬ 
vals  of  10°C  in  the  surface  tanperatrure  decrease. 


For  the  mu  soot  parameters,  results  obtainsd  witn  the  EBM  are  shown 
in  Figure  4  (far  spring  conditions)  end  Figure  S  (for  sunner  conditions) . 
The  soot  effect  mbs  siau luted  by  reducing  the  solar  heating  In  the  30°  to 
70^f  latitude  band  to  5  percent  of  the  Twrwi  value  for  the  season  (that 
is.  s"3  »  0.09).  Results  are  shorn  for  day  10  of  the  siaalations.  Con¬ 
tours  aarfc  intervals  of  2°C  In  the  predicted  tasperature  decrease. 

Ih  cohering  these  aodel  predictions,  the  effect  of  the  finer  spatial 
resolution  of  the  CCM  is  obvious.  The  tasperature  profiles  In  Figure  3 
closely  follow  the  continental  boundarlae.  By  wparlaon.  these  boundaries 
are  blurred  in  the  contours  of  Figures  4  and  5.  Figure  3  also  shams  the 
effects  of  the  large-scale  circulation  In  the  CCM.  For  instance,  cooling 
tends  to  be  less  severe  on  the  western  sides  of  the  continents  than  an  the 
eastern  sides.  This  is  because,  at  aid-northern  latitudes,  the  prevailing 
westerlies  bring  warm  ocean  air  into  the  western  regions  of  the  continents. 
Mo  such  effect  can  occur  in  the  EBM,  *4iich  only  allawe  far  isotropic 
diffusion  of  hast. 

Table  2  cospares  the  mexliaim  tasperature  decreases  predicted  by  the 
two  models.  At  each  position  in  the  table,  the  value  from  the  CCM  is  given 
first,  followed  by  the  corresponding  prediction  of  the  EBM.  This  compar¬ 
ison  show  that  the  EBM  tends  to  underestimate  the  severity  of.  the  tempera¬ 
ture-  decrease,  as  compered  to  the  CMM.  This  may  be  attributed  to  two 
effects.  First,  as  already  noted  in  the  discussion  of  the  normal  seasonal 
cycle,  the  response  of  the  EBM  to  changes,  in  solar  hasting  tends  to  be  too 
sluggish.  Second,  as  a  result  of  the  coarser  spatial  resolution  of  the 
EBM.  local  teoperature  extremes  are  smoothed  out  cr«r  larger  areas  and 
decreased  in  amplitude.  This  comparison  indicates  that  it  would  have  been 
better  to  have  included  more  spatial  modes  in  the  EBM  solution. 


Table  2.  Maximus  temperature  changes  for  the  NCAR  scenario. 


OCM/EBM 

•— 

Western  Hemisphere 

Eastern  Hemisphere 

Spring 

>20°C/13°C 

”30°C/24°C 

Sumner 

>40°C/21°C 

>40cC/3S°C 

4.  Temperature  changes  predicted  by  the  EBM,  following  a  springtime 
•ch)  injection  of  a  uniform  soot  layer  in  a  3(P  to  70°N  band, 
i  to  top  panel  of  Figure  3. _ 


June)  Injection  oj 
Coepare  to  bo  t  ten 


In  spite  of  these  quantitative  differences,  the  EBM  does  reproduce 
nasty  features  of  the  CCM  predictions.  Both  model*  Indicate  that  the 
largest  temperature  drops  are  confined  to  the  deep  continental  Interiors. 

In  both  cases.  Thble  2  shews  that  ths  cooling  la  sere  severe  far  susner 
nuts  and  for  the  Eastern  Hemisphere.  He  conclude  that,  although  the  H3M 
will  not  yield  quantitatively  accurate  predictions,  the  major  qualitative 
features  are  In  agreement  with  GCM  results.  (liven  the  present  (and  near- 
future)  uncertainties  In  Input  data  cn  the  soot  loading*  It  Is  doubtful 
that  any  type  of  global  climate  model  Is  able  to  produce  quantitatively 
accurate  predictions. 

Scientists  at  the  Computer  Center  of  the  O.S.S.R.  Academy  of  Sciences 
(CCAS)  used  a  GCM  to  simulate  the  effects  of  a  very  dense  soot  layer, 
combined  with  a  stratospheric  dust  layer  (Reference  23).  The  CCAS  model 
uses  a  horizontal  grid  with  12°  latitude  by  13°  longitude  resolution  plus 
two  vertical  layers.  Sea  surface  temperatures  are  allowed  to  vary  In 
accordance  with  an  energy  balance  equation  far  the  ocean  mixed  layer.  The 
model  uses  amual  mean  solar  heating,  so  there  Is  no  seasonal  cycle.  The 
soot/dust  parameters  employed  In  the  CCAS  simulation  are: 

(1)  statospherlc  dust  layer  reflects  70  percent  of  the  incoming 
direct  sunlight  back  to  space, 

(2)  tropospheric  soot  layer  is  purely  absorbing  with  optical  depth  * 
6.0  at  visible  wavelengths. 

(3)  latitude  of  soot /dust  »  12°  to  90  °N.  and 

(4)  no  effect  an  Infrared  radiation. 

Figure  6  shows  the  temperature  changes,  after  40  days,  as  computed  by 
the  OCAS  model.  Figure  7  shawm  the  corresponding  results  obtained  using 
the  EBM.  with  the  seasonal  cycle  turned  off. 

The  geographical  distributions  of  temperature  changes  shewn  in  Figures 
6  and  7  are  very  similar.  In  the  Eastern  Hemisphere,  the  major  temperature 
changes  extend  much  further  south  than  in  the  NCAR  scenario  (Figures  3  to 
3) .  For  instance,  large  areas  of  Africa  would  experience  temperature  drops 
of  more  than  10^,  according  both  to  the  CCAS  model  and  to  the  EEM.  This 
occurs  because  the  soot/dust  layer  new  extends  much  further  southward, 
where  the  solar  heating  is  normally  greatest.  Cfcantitatively,  the  CCAS 
model  am  the  EBM  are  in  agreement  with  respect  to  the  maximun  temperature 
decrease  for  the  Eastern  Hemisphere.  The  CCAS  simulation  has  a  larger 
temperature  drop  In  the  Western  ten! sphere  than  the  EBM  prediction  (30°C 
vs.  23°C) .  This  tendency  of  the  EEM  to  underestimate  Western  Hemisphere 
temperature  changes  >ms  also  apparent  In  the  comparison  to  the  NCAR  CCM. 

A  comparison  of  GCM  results  has  been  jointly  published  v  the  NCAR  and 
CCAS  groups  (Reference  24) .  Although  the  two  models  show  gro*.  *  similar 
results,  the  assumed  soct/duat  parameters  and  the  simulated  ell—  ites 
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■e  6.  Temperature  changes  predicted  by  OCAS  coaputei*  Model  (GCM  plus 
ocean  mixed  layer) .  at  40  days  after  the  initial  soot /dust 
perturbation. 


f— Moral  vs.  annual  aeon)  are  too  different  to  allow  any  detailed  com- 
periscns.  This  points  up  an  ambiguity  In  comparing  EBM  results  to  GCM 
results— the  degree  to  which  various  GCMe  would  agree  in  their  predictions 
far  a  given  scenario  Is  not  known. 


SQdTlJITY  TO  ASSQED  SOOT  PARAMETERS 

As  remarked  earlier,  the  major  advantage  of  the  EBM  over  more  complex 
models  Is  caaputatlanal  convenience.  This  allows  far  simulation  of  a  wide 
variety  of  cases.  Although  the  predicted  climatic  effects  should  not  be 
considered  quantitatively  accurate,  carparlacns  of  results  from  different 
scenarios  are  useful  In  showing  sensitivity  to  seamed  parameters. 


The  bsssllne  case  Is  derived  largely  from  the  1985  study  by  the  MAS 
(Reference  7) .  This  scenario  assumes  a  zonal  ly  uni  fare  soot  cloud  confined 
between  latitudes  30°  to  70°N.  The  soot  Is  purely  abearblng  to  Incoming 
sunlight  and  has  no  effect  an  outgoing  infrared  radiation.  Initially,  the 
cloud  has  an  optical  depth  of  3.0  to  sunlight,  one-half  of  the  soot  Is 
placed  in  the  lower  troposphere  (below  5  km) .  there  normal  removal  pro¬ 
cesses  lssd  to  a  half-life  of  3  days  (a- folding  time  of  4.33  days) .  The 
remaining  soot  is  placed  in  the  S  to  10  tarn  layers,  there  the  ha-'f-lius  ror- 
removal  is  30  days  (e- folding  tin  of  43.3  days).  Of  course,  in:  the  one- 
level  EBM,  the  only  effect  of  the  vertical  distribution  of  soot  is  an  the 
assumed  removal  rates.  Symbolically,  the  removal  rates  are  reflected  In 
the  time-dependence  of  the  cloud  optical  depth  at  visible  wavelength: 

-t/4.33  -t/43.3 

Tvi,  <*>  "  •  1 V  *  V  >  171 

where  time  la  now  measured  In  days  since  the  initial  soot  injection  and 

W  "  3-°-  \  -  *o  "  1/2 '  (a) 

.  far  the  baseline  case,  the  soot  Is  assumed  tu  be  injected  during  the  summer 
season  (30  June) . 

Figure  8  shows  Western  Hemisphere  maps  of  the  Induced  surface 
temperature  change  (with  respect  to  normal  seasonal  temperatures),  at  10, 
20.  40.  and  80  days  after  the  soot  injection.  The  maximum  effect  occurs 
near  day  20  of  the  simulation,  **en  surface  temperatures  drop  by  as  much  as 
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23%.  (ftar  sensitivity  studies.  It  is  sufficient  to  consider  only  the 
Hsstsm  Hemisphere,  sines  tbs  other  hemisphere  behaves  in  a  qualitatively 
similar  manner.)  .  .  ^ 

Man  considering  the  effects  of  a  climatic  anomaly  on  agriculture,  the 
degree  day  is  a  cceoonly  need  metric  (Reference  25,  Chapter  7).  Here,  we 
define  degree  days  as  the  difference  between  daily  tam'eiiatnrew  in  the 
perturbed  vs.  control  cliamtes,  summed  over  a  specified  length  of  time. 
Figure  9  displays  a  degree  day  map  for  the  baseline  case,  with  the  sue- 
matlcn  extending  over  a  1-year  period  after  the  initial  aoot  injection.  An 
interesting  aspect  of  this  asp  is  that  the  degree-day  oontours  do  not 
follow  coastlines  nearly  as  closely  as  do  the  temperature  contours.  This 
is  because.  Mills  the  oceans  cool  more  slowly  than  continental  areas,  the 
oceans  also  recover  more  slowly  as  the  soot  is  removed.  TO  a  certain 
extant,  this  decreases  the  effect  of  the  greater  heat  capacity  of  the 
oceans,  in  terms  of  a  time-integrated  quantity  like  degree  days.  TO  the 
extent  that  biological  damage  depends  on  the  duration,  as  well  as  severity 
of  cooling,  it  appears  that  the  consequences  for  serine  biota  are  oempar- 
able  to  those  in  continental  environments,  particularly  since  aarine  life 
forma  tend  to  be  learn  tolerant  to  temperature  changes  (Reference  8) . 

The  sensitivity  tests  Involved  various  modifications  of  the  baseline 
scenario.  These  modifications  and  the  predicted  effects  are  suanarizad  in 
Table  3.  The  remainder  of  this  section  provides  brief  descriptions  and 
cooments  an  these  simulations. 


Seasonal  Dependence 

The  cooling  effect  is  obviously  largest  if  soot  injection  takes  place 
during  the  sunner.  since  that  is  Mien  the  solar  heating  is  normally 
greatest.  With  the  exception  of  the  NCAR/CCM  simulation  (Figure  3) ,  little 
has  been  published  regarding  the  consequences  of  soot  injections  in  other 
seasons.  As  Table  3  shows,  there  is  a  strong  seasonal  dependence,  both  in 
the  maxliwim  cooling  and  in  degree  days.  Perhaps  more  surprising,  the 
effects  are  very  different  for  spring  ve.  autunn  scenarios,  particularly 
in  the  degree  day  column.  This  is  because  the  spring  scenario  occurs 
during  the  phase  of  the  seasonal  cycle  when  solar  heating  and  surface 
temperatures  are  both  normally  increasing,  in  the  northern  Hemisphere,  at 
the  nmxlwum  rate.  The  autumn  scenario  occurs  Mien  surface  temperatures  are 
normally  decreasing,  at  a  rets  controlled,  to  a  large  extant,  by  the  finite 
heat  capacity  rather  than  by  solar  energy  effects. 


Soot  Removal  Rates 

The  baseline  case  uses  the  assumed  altitude  distribution  of  the  soot 
to  estimate  the  removal  rates,  based  on  normal  atmospheric  cleansing 
processes.  Cases  5  to  7  examine  the  dependence  of  the  calculated  climate 


Trial*  3:  Susury  of  Sensitivity  Tuts 


Cae*  0 

Modification  of  Baseline 

Max.  Cooling 
(°C) 

Degree 

-Days* 

1. 

Basellne/Susner 

22.8 

2388 

2. 

Baeellne/Sprlng 

17.0 

2285 

3. 

Baseline/Winter 

6.9 

1092 

4. 

Basel lne/Autuan 

12.5 

1188 

3. 

* I-  XtJ  -  0 

16.3 

384 

6. 

^  -  2/3.  Xg  «  1/3 

19.6 

1901 

7. 

*L  ■  0'  *0  *  1 

28.8 

3322 

3. 

Narrow  latitude  band 

20.5 

2147 

9. 

Two  latitude  bends 

14.1 

1733 

10. 

Outgoing  infrared  reduced  50* 

14.0 

1322 

a 

Summed  over  one  year  period  after  Initial  soot  Injection. 
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perturbation  on  the  soot  removal  rata.  In  Case  5.  tbs  soot  la  summed 
occur  only  In  the  loner  troposphere.  v*er»  rapid  removal  occurs.  The 
climatic  perturbation  Is  then  less  severe  and  of  very  short  duration,  as 
indicated  by  the  factor  of  six  reduction  In  degree  days.  In  Case  6.  the 
soot  is  partitioned  between  the  loner  and  upper  troposphere  In  the  ratio 
which  represents  a  constant  mixing  ratio  in  the  0  to  10  km  altitude  range, 
m  case  7.  all  of  the  soot  occurs  In  the  upper  troposphere,  there  the 
lifetime  Is  very  long. 

A  caution  regarding  Interpretation  of  these  results  Is  in  order.  Since  the 
EBM  is  a  ana-level  model,  it  cannot  treat  the  effects  of  the  scot  altitude 
distribution  on  radiative  processes.  Ramaawry  and  Klehl  (Reference  26) 
used  a  one-dimensional  (vertical)  radiative  convective  model  to  examine 
atmospheric  cooling  rates  for  different  soot  altitude  profiles.  They  found 
that  the  surface  cooling  rate  is  cut  In  half  If  a  constant  mixing  ratio 
with  a  scale  height  of  3  km  is  assumed.  This  occurs  because  much  of  the 
solar  radiation  Is  absorbed  user  the  planetary  boundary  layer  which  Inter¬ 
acts  strongly  with  the  surface.  Such  effects  are  neglected  In  the  HEM. 
Furthermore,  the  assumption  of  a  fixed  height  distribution  (modified  only 
by  removal)  Is  very  artificial.  In  reality,  the  radiatlvely  heated  soot 
cloud  will  tend  to  rise  over  an  extended  period  of  time,  as  discussed  in 
Section  3. 


Horizontal  Soot  Distribution 

In  the  first  few  days  after  a  nuclear  war.  the  scot  would  be  confined 
to  isolated  clouds  or  plumes,  over  the  fire  sites.  If  the  soot  remains  in 
the  atmosphere  for  several  days,  normal  circulation  plus  winds  created  by 
horizontal  thermal  gradients  vwuld  tend  to  move  the  soot  way  from  the  fire 
sites.  The  geographical  distribution  of  the  soot  cloud  would  then  depend 
on  the  distribution  of  fire  sites,  on  normal  circulation  patterns,  and  on 
the  anomalous  circulations  created  by  locally  uneven  heating  and  cooling. 

It  is  by  no  means  obvious  that  the  result  of  transport  would  be  a  cloud 
layer  v*iich  is  approximately  uniform  in  a  broad  latitude  band. 

Because  climate  models  are  barely  able  to  resolve  the  synoptic  scale, 
climate  simulations  have  assumed  that  the  scot  cloud  is  uniform  between 
specified  boundaries.  The  simplest  case  (disregarding  the  totally  unreal¬ 
istic  case  of  a  global  cloud)  is  a  zonal iy  uniform  cloud  between  specified 
northern  latitudes.  This  assumption  may  be  partially  justified  by  the 
similar  latitudes  of  the  primary  target  areas  in  both  the  Eastern  and 
West  era  Hemispheres .  and  by  the  highly  zonal  character  of  normal  circula¬ 
tion  patterns  at  the  target  latitudes.  For  instance,  the  baseline  case 
assumes  that  the  soot  lies  In  a  band  between  latitudes  30°  to  70°N,  as  did 
the  NCAR  study  ( Reference  20 ) . 
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1b  tart  the  sensitivity  of  the  Billeted  ollaete  to  this  assumption. 

8  and  9  used  zonal  cloud  bands  with  different  boundaries.  For 
8.  the  cloud  bend  wee  confined  between  40°  and  80%.  Since  this 
only  one  half  the  cloud  extant  In  the  baseline  case,  the  optical  depth  wee 
doubled  (l.'a. .  x_  a  6.0)  so  that  the  asaueed  total  soot  ease  was  the  sane 
In  both  eases.  JBnxilcated  In  Table  3.  this  results  In  a  very  slight 
reduction  In  the  severity  of  the  eliaate  effect,  itots  that.  In  Case  8*  the 
optical  depth  under  the  cloud  Is,  at  any  given  tine,  twice  as  large  as  in 
the  baseline  case.  As  a  result,  the  optical  depth  is  greater  then  unity 
for  48  days  in  Cues  8,  as  opprjsal  to  18  days  in  the  base  line  cess,  fills 
longer  duration  is.  however,  completely  compensated  by  the  fact  fiat  the 
center  of  the  cloud  band  is  now  such  closer  to  file  unshaded  latitudes,  so 
horizontal  diffusion  of  beet  fra*  these  latitudes  is  sore  Important. 


9  takes  s  different  approach  to  the  Issue  of  eoot  cloud  location. 
Bis  soot  was  confined  to  two  bends,  between  20°  to  40%  and  between  60°  and 
80%.  The  solar  heating  In  the  gap  (40°  to  60%)  wee  left  at  normal  values 
to  elanlate  file  effect  of  a  broken  cloud  layer.  Tha  sexism  temperature 
effect  fear  this  case  Is  40  percent  smaller  then  for  the  besellne 
while  the  degree  dey  effect  is  reduced  by  30  percent  (ass  Thble  3) . 


The  geographical  distribution  of  the  climate  ancaaly  far  this  case  is 
very  significant.  Figure  10  shews  a  temperature  difference  map  (Case  9 
minus  control)  an  day  30.  which  corresponds  to  the  greatest  cooling;  the 
degree  day  distribution,  sunned  over  1  year,  is  shown  In  Figure-  11.  Note' 
that,  in  both  figures,  tha  greatest  excursions  occur  in  the  bend  covered  by 
the  more  northern  soot  cloud.  This  is  surprising  in  that  the  larger  abso¬ 
lute  reduction  In  solar  heating  occurs  In  tha  more  southern  soot  band, 
where  the  noised  solar  flux  Is  larger.  As  in  Case  3.  horizontal  transport 
of  heat  plays  a  very  Important  role  In  these  simulations  with  narrow  cloud 
tmnds.  The  southern  soot  bend  is  hasted  an  both  sides  by  horizontal  heat 
transport  from  unshaded  latitudes.  The  northern  soot  bend  extends  nearly 
to  the  pole.  As  a  result,  the  unshaded  region  northward  of  this  band  Is 
email  In  ares  and  receives  little  solar  energy,  sc  the  region  In  the  soot 
band  can  only  be  heated  by  horizontal  transport  from  the  unshaded  area  to 
the  south.  This  simulation  points  out  the  importance  of  heat  transport 
from  clear  areas  into  soot  regions  vixen  the  soot  cloud  is  limited  in 
extent.  For  reasons  given  above ,  the  southern  soot  band  say  be  more 
representative  of  this  effect.  In  this  region,  the  smxlmum  temperature 
Is  only  "  10%,  l.a. ,  lass  than  ore-half  of  the  temperature 
in  the  baseline  case. 


As  argued  earlier,  there  are  no  compelling  reasons  to  expect  that  the 
soot  cloud  would  expand  into  a  broad  and  uniform  zonal  band.  Tracer  calcu¬ 
lations.  using  normal  atmospheric  circulation  as  Input,  indicate  that  the 
soot  cloud  will  spread  exit  frem  the  target  areas,  but  this  spreading  would 
be  far  from  uniform  (Reference  27) .  This  calculation  probably 
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Figure  id.  Temperature  differences  (perturbed  minus  control)  at  day  30  in 
the  Case  9  simulation. 
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Figure  11.  Degree  day  map  for  the  Case  9  simulation, 
period  following  the  initial  smoke  injection. 


sunned  over  a  1-year 


overestimates  the  degree  of  uniformity,  doe  to  limited  epetial  resolution 
of  the  input  winds  nl  due  to  nsnsricel  diffusion  In  the  transport 
algorithm. 

Zt  would  be  of  Interest  to  simulate  the  effect  of  narrower  soot  bands 
than  considered  In  Cases  8  and  9.  as  well  as  soot  clouds  having  longitu¬ 
dinal  'variations,  in  addition  to  latitudinal  structure.  Longitudinal 
variations  would  allow  hast  transport  to  occur  in  both  of  the  horizontal 
dimensions,  this  was  not  passible  with  the  current  version  of  the  EEM 
because:  (a)  the  limited  number  of  modes  retained  in  the  spherical 
harmonic  expansion  [Equations  (4)  to  (6)1  sets  the  resolution  at  about  15° 
of  latitude:  and  (b)  the  current  implementation  of  the  EEM  assumes  a 
zonally  uniform  soot  cloud.  However,  neither  of  these  limitations  are 
Inherent  In  the  basic  model,  so  these  Issues  could  be  examined  In  future 
Investigations . 


Outgoing  infrared  Radiation 

The  climate  simulations  assume  that  the  soot  cloud  is  transparent  to 
Infrared  radiation.  This  is  a  critical  assumption  since  absorption  of 
incoming  sunlight,  coupled  with  the  transparency  to  outgoing  thermal 
energy,  leads  directly  to  the  predicted  cooling.  While  It  is  true  that 
subnicron  particles  of  pure  soot  absorb  mainly  at  visible  (or  shorter) 
wavelengths,  the  actual  cloud  resulting  from  fires  would  not  consist  of 
pure  soot.  Large  amounts  of  water  vapor,  carbon  monoxide,  and  nitrogen 
oxides  would  be  entrained  in  the  fire  pluses,  along  with  a  vast  inventory 
of  hydrocarbons  and  toxic  chemicals  produced  by  incomplete  combustion 
(Reference  28) .  Many  of  these  gases  are  strong  absorbers  in  the  infrared. 

To  test  the  effect  of  increased  Infrared  opacity,  we  reduced  the 
outgoing  Infrared  fluxes  In  the  fflM,  within  the  soot  band.  The  maximum 
reduction  was  30  percent,  with  a  return  to  normal  fluxes  as  the  soot  was 
removed.  As  expected,  the  severity  of  the  climate  modification  was 
significantly  reduced  (see  Table  3).  Ohfortunately,  a  realistic  estimate 
of  the  infrared  opacity  will  require  detailed  consideration  of  atmospheric 
chemistry,  which  is  beyond  the  scope  of  this  report.  The  results  of  Case 
10  show  that  the  Nuclear  Winter  problem  is  sensitive  to  Infrared  opacity, 
but  the  predictions  must  not  be  interpreted  as  a  quantitative  estimate. 
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Tho  EBM  simulations  rVurrlbad  In  Section  2  doeonetrate  the  appli¬ 
cability  of  this  type  of  nodal  to  the  soot-cooling  problem.  Wills  the 
results  should  not  be  interpreted  as  quantitative  predictions,  they  are 
useful  in  delineating  sensitivity  of  the  predicted  cl  lasts  modification  to 
uncertainties  in  the  input  parameters.  The  more  ccaqalex  GCMs  share  these 
uncertainties  in  that  radiation,  scavenging  processes .  tropospheric  che¬ 
mistry,  and  sub-grld-acale  transport  are  either  heavily  parameterized  or 
totally  Ignored.  Such  models  are  too  cumbersome  to  allot*  the  type  of 
sensitivity  testing  described  in  Section  2.  Thus,  the  flexibility  of  the 
EBM  to  simulate  many  cases  complements  the  more  powerful  predictive 
capabilities  of  the  GCMs. 

The  important  new  implication  to  be  drawn  from  the  sensitivity  tests 
is  that  the  severity  of  the  climate  modification  depends  critically  on  the 
lateral  extent  and  uniformity  of  the  soot  cloud.  If  the  cloud  is  broken  or 
confined  to  narrow  strips,  hast  frcr  nearby  unclouded  regions  will  be 
transported  into  the  clouded  areas,  thereby  limiting  the  surface  tempera¬ 
ture  decrease,  there  has  been  speculation  that  a  very  nonuniform  soot 
layer  would  have  a  reduced  climatic  iapact  because-  the  soot  particles  would 
be  less  effective  in  absorbing  sunlight  (l.e. ,  optical  depths  would  vary 
between  very  large  and  very  small  values) .  While  this  reasoning  is 
correct.  Case  9  shows  that  the  climatic  impact  is  reduced  even  if  the 
broken  soot  cloud  is  distributed  such  that  the  total  sunlight  absorbed  is 
equal  to  that  absorbed  in  the  baseline  case  (l.e.  a  broad  and  uniform 
cloud) . 

This  case  points  to  the  importance  of  properly  accounting  for  the 
lateral  transport  and  spreading  of  the  soot  cloud  away  from  localized  fire 
sites.  While  the  basic  physics  of  atmospheric  transport  is  well  under¬ 
stood.  the  ability  to  model  this  process  is  limited  because  such  modeling 
requires  both  high  resolution  (down  to  the  plane  scale)  and  large  domain 
(up  to  the  planetary  seals) .  GCMS  will  definitely  not  be  able  to  properly 
handle  this  problem,  so  appropriate  parameter izatiens  of  sub-grid-scale 
soot-cloud  structure  will  have  to  be  developed.  Simpler  models,  such  as 
the  EBM.  can  be  useful  in  developing  and  testing  various  schemes  for 
simulating  the  effects  of  nonuni  form  or  broken  soot  clouds. 

The  EBM  wes  also  used  to  demonstrate  sensitivity  of  the  climata  effect 
to:  the  assumed  season,  the  soot  removal  rate,  and  the  effect  of  the 
atmospheric  loading  an  infrared  radiation.  In  these  cases,  the  model 
results  confirm  intuitive  expectations  that  the  severity  of  the  climate 
modification  is  very  dependent  on  these  factors. 
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As  s  final  note,  ws  offer  a  warning  regarding  interpretation  of  aodel 
predictions  for  the  cliaate  aodification  scenario.  Men  dealing  with  very 
coepiax  at. dels,  such  as  GOb,  it  is  taspting  to  accept  the  predictions 
without  scrutinizing  **sther  the  basic  processes  are  properly  handled  by 
the  models:  the  trnaendous  amount  of  detail  in  the  sodels  can  give  a  false 
sense  of  security.  An  e»nple  will  be  provided  to  illustrate  this  danger. 
based  on  the  recant  GO!  simulation.  carried  out  at  I cm  Alans  (Reference 
29) .  Ibis  simulation  represents  the  current  state-of-the-art  in  global 
effects  Bodeling. 


sEXJHwrrxNG  or  the  socrr  laser 

The  major  improvement  of  the  Los  Alanos  simulation  over  other  models 
Is  that  thn  soot  Is  Interactive  (i.e. .  directly  participates  in  the  energy 
balance  and  dynamics  of  the  simulated  atmosphere) .  The  most  significant 
new  feature  to  emerge  from  the  Los  Alamos  simulation  is  that,  because  the 
soot  absorbs  sunlight,  it  heats  and  rises  into  the  stratosphere.  This 
self-lofting  is  critically  important  because  the  lifetime  of  the  soot  is 
much  longer  once  the  soot  rises  above  the  precipitating  layers  of  the 
atmosphere. 

There  is  little  question  that  there  is  a  potential  for  self-lofting  if 
the  soot  layer  survives  initial  scavenging.  However,  it  is  doubtful  that  a 
GCM  can  directly  follow  this  process.  Figure  12  provides  a  schematic 
illustration  of  how  the  lofting  process  will  likely  proceed..  The  top  panel 
shows  an  initial  configuration  consisting  of  a  horizontally  uniform  soot 
layer.  Since  the  optical  depth  is  greater  than  unity,  much  of  the  3oot  is 
shielded  from  direct  sunlight.  The  middle  panel  shows  a  perturbation  of 
the  uniform  cl  cud.  consisting  of  a  local  upward  intrusion.  Such  a  pertur¬ 
bation  could,  for  instance,  be  created  by  local  convective  processes.  This 
intrusion  will  be  heated  more  rapidly  than  the  remaining  soot  layer  because 
it  provides  a  larger  surface  area  for  intercepting  sunlight.  As  a  result, 
the  intrusion  will  rise  faster  than  the  uniform  layer,  as  illustrated  in 
the  bottom  panel. 

This  instability  of  the  uniform  soot  layer  is  an  example  of  multi¬ 
component  (or  "doubla-dif fusive” )  convection  (Reference  30).  Such  pro¬ 
cesses  have  been  widely  studied  in  the  laboratory  and  in  oceanography, 
where  the  intrusiare  are  referred  to  as  "salt  fingers".  Figure  13  shows  a 
laboratory  example  of  up»mrd-propagatlng,  double-diffusive  convection. 

In  contrast  to  this  highly  structured  physical  process .  GCMs  are 
restricted  to  the  .■ery  uniform  clcud  layer  shown  in  the  top  panel  of  Figure 
12.  The  Lea  Alamos  GCM,  which  is  based  an  the  NCAR/CCM,  has  horizontal 
resolution  of  roughly  500  km  and  vertical  resolution  of  the  order  of  1  km 
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Figure  12.  Development  of  a  double-diffusive  instability  In  a  uniform  eoot 
cloud.  Tog  panel:  initial  conditions,  consisting  of  a  horizontally 
uniform,  optically  thick  cloud.  Middle  panel:  perturbation  of  the 
initial  condition.  Bottom  panel:  growth  of  the  perturbation  due  to 
enhanced  solar  heating. 


Figure  13.  A  laboratory  erample  of  double-diffusive  convection  produced  by 
placing  hot  salty  water  over  cool  fresh  water.  The  fresh  water  is 
marked  by  flourescent  dye.  Fran  Reference  30. 


at  the  cloud  level.  In  other  words,  the  beslc  aodel  call  is  a  very  thin 
plats  (aspect  ratio  of  900:1) .  Sjch  a  aodel  cannot  directly  follow  the 
self-loftlng  process,  tdilch  will  have  a  horizontal /vertical  aspect  ratio 
«1  (as  In  Figure  13).  In  fact,  the  doubla-dlffuelve  Instability  aalces  1** 
uni  Iks  ly  that  a  broad  and  horizontally  uniform  cl  nod  will  ever  farm. 
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BXHDMICAL  Httgaus  or  IB  CLIMATE  M3CBL 


Ha  eliti  aodtl  defined  by  Equation  (1)  has  both  tins-  and  apses 
dcpmdmct.  Has  Mathematical  propertlsa  of  this  aodel  concern  the 
•-folding  tins  for  response  to  a  time  dependant  perturbation,  and  the 
characteristic  length  for  propagation  of  a  localized  perturbation.  In 
principle,  these  properties  could  be  derived  by  applying  a  formal  Green's 
function  analysis.  However,  tills  is  unnecessarily  cceplesc;  mb  take  a  less 
rigorous  approach  based  an  escaninatlan  of  ths  specific  terse  governing 
tine-  and  space- propagation  of  a  perturbation. 

1b  derive  the  e-folding  tine,  we  consider  a  spatially  averaged  version 
of  Equation  (1).  The  diffusion  tom  then  vanishes,  leaving 


QeS(t)  -  CA  ♦ 


(A.l) 


Me  consider  S(t)  to  be  the  source  term  far  driving  a  step-function 
perturbation  of  this  equation: 


tSQ  for  t  <  0 
3^  for  t  >-  0 


(A.2) 


Assuming  that  ths  modal  Is  in  equilifaritm  for  t  <_0.  at  a  temperature 


To  *  (<JaSo  "  A1/B- 

the  solution  to  Equation  (A.l)  far  t  >  0  is 


(A.3) 


T(t)  »  Tj  +  (T0  -  Tj)' 


Tt  -  (QaSj  -  A)/B 


-Bt/C 


Frcm  Equation  (A.4) .  we  see  that  T(t)  ->  T..  as  t  ->  »,  and  that  the 
•-folding  time  for  the  temperature  response  Is 


1b  derive  the  characteristic  laagth  scale,  mb  rewrite  Equation  (1)  in 
the  farm 


I9T  B  OaS  -  A  012 

- -  -  ( - 1)  + - 7  T 

Ijt  C  ST  C  T 


(A.7) 


ttoe re  time-  and  space-  dependence  of  the  appropriate  variables  is  still 
implied.  Examining  this  equation,  it  is  apparent  that  the  left  side  is  the 
inverse  of  the  e-folding  time  for  teaperature.  The  first  term  an  the  right 
side  involves  cnly  radiation  quantities  (solar  and  infrared),  while  the 
second  ter j  represents  horizontal  energy  transport.  The  teaperature  time 
scale  is  thus  the  harmonic  sum  of  the  radiation  and  transport  time  scales: 


Ttenperature  Tradiaticn  1  transport 


(A.  8) 


where  the  time  scales  are  evaluated  by  term-by-term  comparison  of  Equations 
(A.7)  and  (A.8) .  He  now  consider  the  consequences  of  a  localized 
perturbation  in  temperature.  EUe  to  the  diffusion  term,  the  perturbation 


will  propagate  outward  at  a  rate  set  by  t, 


transport 


.  As  it  propagates,  it 


will  be  damped  by  radiative  processes.-  at  a  rate  set  by  t  .  .  ^  \ 

characteristic  distance  to  vtfilch  the  disturbance  will  propagate  is  obtained 
by 


radiation 


D  1  B  OaS  -  A 

- 7  T*  »  -  ( - 1) 

C  T  C  3T 


(A. 10) 


The  characteristic  length,  l.  at  the  diffusion  operator  is  given  by 


1  1  B  OaS  -  A 

—  -  -  7  2T  -  -  { - 1) 

f  T  D  BfT 


(A.  11) 


Finally,  since  the  equilibrium  teaperature  is  given  by 
QeS  -  A 


(A. 12) 
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